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Abstract 



We calculate the dynamical behaviour of a hole in various spin back- 
grounds in infinite dimensions, where it can be determined exactly. We con- 
sider hypercubic lattices with two different types of spin backgrounds. On 
one hand we study an ensemble of spin configurations with an arbitrary spin 
probability on each sublattice. This model corresponds to a thermal average 
over all spin configurations in the presence of staggered or uniform magnetic 
fields. On the other hand we consider a definite spin state characterized by the 
angle between the spins on different sublattices, i.e a classical spin system in 
an external magnetic field. When spin fluctuations are considered, this model 
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describes the physics of unpaired particles in strong coupling superconductors. 



Typeset using REVTgX 



I. INTRODUCTION 



It is well known that the discovery of high Tc superconductors has triggered an extensive 
study of highly correlated systems. The Hubbard and t — J models have been prototype 
models towards the understanding of most of the features of those materials. For low doping 
(near half filling) the system consists of dilute mobile holes in a spin background. 

A lot of effort has been devoted to the understanding of the dynamics of holes in spin 
backgrounds. Brinkman and Rice (BR) || considered several configurations of the spin 
background (ferromagnetic (FM), Neel and random) and studied the density of hole-states 
and dc conductivity. Their calculation was based on Nagaoka's expansion of expectation 
values in terms of background-conserving hole paths |§. Of these, they only considered 
retraceable paths, i.e., no loops were taken into account. They obtained very accurate 
results for single-particle Green's functions of a hole in the Neel background, being exact 
in one dimension (where all paths are retraceable). This retraceable paths approximation 
(rpa) was used to study many other quantities such as dynamical conductivity [|J, electrical 
resistivity, thermal conductivity, thermopower and specific heat [[J. 

Using different approximations, many other analytical approaches have been reported 
for the ground state and excited properties of a single hole |||7|,[J . There are also numerous 
studies using exact diagonalization techniques in low dimensions P, |TU| , |rT| , |T2| . These results 
show a well defined quasiparticle peak for J ~ t whereas for small enough J ft an incoherent 
spectrum carries most of the spectral weight. 

Most of the studies have been concentrated on the dynamics of a hole in an antifer- 
romagnetic (AF) background. The case of the polarized t — J model has been studied in 
connection with strong coupling superconductors [|H|[13|. The negative U Hubbard model 
with n particles is equivalent to a positive-?/ case with one particle per site and a net mag- 



netization given by S z = ( - 1 ~ n - ) ||14|| . So the dynamics of a single unpaired particle moving in 
a background of strongly bounded paired particles (strong negative U limit) is described by 
the t — J model. 
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In the present work we study the single-particle Green's functions in certain backgrounds, 
extending the results of Metzner et al. |IJ to generalized spin backgrounds and in particular 
to one that describes the physics of holes in strong coupling superconductors. 

The calculations are performed for infinite dimensions where the results are exact. The 



limit of high lattice dimensions, d — ■> oo, has been used to study correlated fermions [T3 



and helped to clarify the validity of several approximations and construct new ones [16[ . 

A self-consistent approximation for finite dimensions can be performed in a similar way 
as in Ref . |l| . More realistic calculations applicable to high Tc superconductors, for example, 
should also include interactions between holes and spin fluctuations. We neglect spin fluc- 
tuations since they disappear in infinite dimensions, and concentrate on the corrections due 
to the inclusion of loops to the BR retraceable path approximation. This latter approach 
doesn't distinguish between different spin backgrounds since paths without loops are always 
background conserving. While in the Neel background the BR approximation is correct 
up to order l/d A , where d is the dimension, the contribution of loops becomes important 
whenever there are clusters of aligned spins. In particular, in a FM background, any hole 
path leaves the background unchanged, leading to Nagaoka's theorem ||. 

The paper is organized as follows: in Sec. II we calculate the dynamics of a hole in an 
ensemble of spin configurations in a hypercubic lattice, considering arbitrary spin averages 
or probabilities in each sublattice; in Sec. Ill we consider a definite spin state characterized 
by the angle between the spins on different sublattices and we obtain expressions for the 
local and k— dependent propagator. We summarize in Sec. IV. 

II. ONE HOLE IN A GENERALIZED SPIN BACKGROUND 

We consider the t — J model for infinite dimensions with J = 0, which with standard 
notation reads: 

H = -t Yl (1 - n^ a )c\ a c jcr (l - n^ a ) (1) 

<ij>,a 
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To keep the average kinetic energy finite in d — > oo, one must scale the hopping amplitude 

t as 

t* 



t* fixed, where Z is the number of nearest neighbours (Z = 2d in a hypercubic lattice) ||15|| . 
The Green's function of a hole in an arbitrary spin background S is: 



Gfj( z ) — X^yV 



z 
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G s ijc {z) = (d—^c^s (3) 

where (• ■ -)s — J2s' w s'(S'\ ■ • • \S') and Ws> is a normalized distribution of spin configurations 
S'. The spins of different sites are statistically independent. We consider a hypercubic lattice 
and characterize the spin ensembles by asigning a probability for spin o in the sublattice 
X: pxa] X = (A,B) and a = (t,l)- The random case considered in Ref. corresponds 
to pAa = Vb<j and interpolates between the FM {j>Aa = 1) and the unpolarized random 
(PAa — 1/2). We generalize the calculation to interpolate also between the random and the 
Neel background {pAo = Pb-c — 1)- The unpolarized random case corresponds to a thermal 
average over all spin configurations which contribute with equal weight. The other polarized 
cases correspond to thermal averages in presence of uniform or staggered magnetic fields. 

We follow closely Ref. for the notations. The calculations are based principally on 
Nagaoka's expansion for the Green's functions in powers of t/z 0. For our hypercubic 
lattice we have to distinguish two local Green's functions: 

Gaa(z) = — 1 and G BB (z) = — 1 (4) 

z[1-Sa{z)\ z[1-Sb{z)\ 

where Sx {X = A, B) is given by the sum over all graphs for which the hole returns to the 

starting point only once. Due to the simple topology of the loop trees, Sx can be written 

as a sum over loops with dressed vertices Cx- This dressed vertex is given by the bare one 

plus all possible Sx insertions, so we have: 

CxA'A = px, T,[Sx(z)} 1 = Vx ° A (5) 



and 



a n=l 

oo 



* s 2n 



<t n=l 



2 

l* \ 2n 



where u n is the number of self-avoiding return paths of length n. 
Defining a generating function for u n 



n=2 



equation (0) can be written as 



Cact 



Sb(z) = E 



M[^C Aa C B(T } - 1 



From Eqs. (|]) and (||) and defining G° as the free particle Green's function 
its inverse, the following relations hold in the FM limit {pa\ = Vb\ —* 0) 



M (0 = 1 + 0(£ 2 ) for small £ M(tG°) = zG c 



This implies 



M(0 = 1^(^/0 



Using this relation in a similar way as in |T| we find: 



1 +PAlPAl(zG A A ~ 1) = ^VA-a\jVAaVBaGAAGBBG\ 



PAaPBaG aaGbb 



1 + Vb\Vb\(zG B b - 1) = ^Vb-cj^VAvVbuG aaG bbG\ 



VAoVBaG aaGbb 



A more compact expression can be obtained in the symmetric case where pAa 
This implies Gaa = Gbb = G and Eqs. ( [TT]) simplify to 

'1 + 7 2 (^G- 1)" 



7 G = G° 



7 G 
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with 7 = ^/pAaPA-cr ranging from Neel to uniform random backgrounds (0 < 7 < 1/2). The 
density of states D(u) = —^ImG(uj + i0 + ) for this symmetric case coincides exactly with 
the density of states of the tilted configuration with the correspondence 7 = a for a < 1/2 
(see Eq. fl2|) and Fig. 2). 

As for \u\ > t* G°(oj) ~ 1/w, Eq. (0) implies ~ 1/w and 7 G ~ G°(uj/j). 

Using the density of states of a free particle at large \u\ we have: 

D{u) ~ ^^ e - 2 /(27 2 ^ 2 ) (13) 
V27T7t* 



The density of states has exponential tails for large \u\ whenever 7 7^ 0. These come from 
the presence of FM clusters in the spin configuration. 

In the Neel limit (pai — > =>- 7 — > 0), to first order in pai, Eq. ( pi]) leads to 



zG - 1 = t* 2 G 2 =► G = V 9 (14) 

2t* 2 



This coincides with the results of the rpa 0. The corresponding density of states 

1 

27rf* 



D M = 7^2^ 4t * 2 -^; H<2t* (15) 



has a semielliptic shape with band edges at ±2t* with a square root singularity. 

We also calculated the self-energy for this symmetric case. It is /c-independent for d — > 00 
and from symmetry E^ CT = E^ ff . It is not necessary to use the Nagaoka expansion for the 
off-diagonal propagator in the site representation as we did for the calculation of the local 
propagator jl]. Instead, we use the Dyson equation together with the locality of the self 
energy. 

In a matrix representation for the fc-dependent Green's function in the reduced Brillouin 
zone ( 6k = — 2i(cos &! + ••• + cos k d . < 0) we have 



( z £ e \ 

I 6fc 2 — Tj Bcr J 



(16) 



This leads to 
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G k(T (z) - — 



Z — S^Bcr — Cfe 
— € k Z — Y>Aa 



D 



z - T, Aa ) (z - S Bct ) - e k ^j (^J (z - Y, Aa ) (z - E B(7 ) + e k ^j 



After some algebra and considering the fact that 



(—^—)k 
a + e k ' k 



o D°(e) r°° D°(e 



a + e 



a — e k 



a - e k a + e k 
and that Gxa = (G\«j)k — VxaG we obtain: 



o a — e 



a - e fc 





— ^Bcr 




— ^Aa 


lz 


— S^o- 



G° 
G° 



z - T, Aa )(z - T, Ba ) 



V z - S B(7 
From Eqs. fll2D and (p0| ) we obtain: 



,2 — £U<r)( z — E S(T ; 



PAaG 



Y, Ba = z(l - p Aa ) 



1-7 2 
1-7 2 



(17) 



(18) 



(19) 



(20a) 
(20b) 



(21a) 
(21b) 



These expressions lead to the known results for the Neel limit (7 = 0) and the random 
case (7 = 1/2) (see Ref. |T[). Here ImT.xaiyJ — «0 + ) > \/u and there are no quasiparticles 
in the system. This is expected since this system doesn't have a Fermi surface. 



III. ONE HOLE IN A TILTED SPIN BACKGROUND 

A more compact and straightforward calculation can be done by considering a magnetic 
field acting on the system polarizing the spins in a way sketched in Fig. 1 . As in the 
previous section, we are considering a hypercubic lattice in a pure spin configuration with 
spins a A and Ob lying on the yz plane, forming an angle betwen them. In this way we can 
continuously connect the Neel (<ft = n) to the FM (0 = 0) background. 
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As we mentioned in the Introduction, this model also describes the dynamics of an 
unpaired particle in a strong coupling superconductor. The spin variables up and down 
correspond to an empty and a doubly occupied site, respectively. The external magnetic 
field plays the role of a chemical potential since the total magnetization is given by the total 
number of particles S z = (1 — n)/2. The magnetization in the x — y plane corresponds to 
the superfluid order parameter. 

By moving along a loop, a hole will replace A-spins with S-spins and the overlap between 
two states at any given site is 

« = (T|e^ /2 |T) = cos(^) (22) 

For a 2n-loop, 2n — 2 sites are changed, therefore 

S(z) = ^u 2n C a (-J = — (23) 

where M is given by Eq. (^). 
Here again 



GW =w=m a,ld c{z ^t^w) < 24) 



From Eqs. (|T0|), (||) and (||) we obtain 

I -a 2 



aG(z) = G° 



az + 



(25) 



aG(z 

This expression is equal to Eq. ([12]) for the symmetric random model with the corre- 
spondence 

a = COs(0/2) = 7 = y/pAaPA-a (26) 

for a < 1/2. But Eq. ( p5|) is more general because a can assume any value between and 
1: a = 1 corresponds to the FM case and a = to the Neel case. The hole dynamics for 
a = 1/2 (=^> = 27r/3) happens to coincide exactly with that for the unpolarized random 
background. 



In Fig. 2 we show the density of states corresponding to the one-particle Green function 
of Eq. (53) for several values of <j). G has been calculated self-consistently using the Neel 



Green function (fH| ) as a seed. The only case without tails is the Neel configuration (0 = tt) 
that presents square root singularities at the band edges. The states in the tails of the other 
configurations correspond to low energy (say uj — > — oo) states due to FM clusters. It also 
seems that for all values of a there is an energy u Q such that D{ojq) is independent of a 
(this crossing also happens for the real part of the propagator at another energy). But this 
coincidence is only approximate, being, nevertheless, a curious feature. 

We calculated the full propagator in the tilted background. For the disposition of axes 
chosen, the spin on the A-sites can be written as: 



^/*|T>=cos(?)|T>-sin(?)|l> 



'4 



(27) 



and on the 5-sites by changing by —(f). So the Green function can be written in a matrix 
way in the spin representation 



/ 



Gaao 



G 



cos 



-sm^jcos^ 



-sm^jcos^ 



sm 2 (f) 



G 



BBo 



G{<j> <- -0) 



(28a) 



(28b) 



For this type of configurations there is a generalized Bloch theorem. The elements of the 
symmetry group that leave the configuration \X) and the Hamiltonian invariant are g\ = 
Tie llna ~ z / 2 , where / = J2t=i ^nj *- e - translations coupled with rotations in 7r about the z 
axis. So cj^ = -7= Xa c^e^' 1 is an eigenoperator of G\ with eigenvalue q such that e* q l = 
(ai) l e l]i ' 1 = e^+^Q/ 2 )- 1 an d Q = (n, - ■ - it). This implies that cj^ can couple only to cj c+Q _ CT . 
We define a propagator in the spin representation 

( ((<bM) ((Cktl c k+Qi)) 

G-kacj' — 

\((c Ck+Q||Ck+Q|)) J 

Then 



(29) 
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(30) 



Considering the locality of the self energy we can write the inverse of Eq. (|29| ) as: 



Gkaa> ( Z ) 



z — Cfc — E 



n 



'TJ 



z + e k - E 



\ 



U / 



so 



Gko-cr'(^) — — 



2 + e fc - Ejj S;^ 
Efj. ^ — efc — S| 



\ 
/ 



D = (z-e k - S TT )(2; + e k - E u ) - E n E iT = -(e fc - £+)(e fc - £L) 



where 



This leads to 



Ckcrcr' 



— .EL 



-.-/•' -S n z-E + -E TT 



Using Eqs. (gDj) and §5§ we obtain 
, /r. x 



Gaao 



(31) 



(32) 



(33) 



E± = - S TT ) ± -v/(S u - S TT )2 - 4E n E iT + A(z - E TT )(z - E u ) (34) 



(35) 



z _i (Ett+Eu )] (G o„ G (| ) _i v ^ (g o +g o ) £ iT (G°_-G^) 
y E n (G°_-G°) ^-i(S TT +E u )](G»-G ( |)+iv / W(G«+G° + ) J 

(36) 



where 



W = (E u - E TT ) 2 - 4E Ti E iT + 4(z - E TT )(z - E 



1 1 / 



(37) 



and we have defined G°(E±) = G±. 



Eq. ( p^) implies S^(z) = Ej,f(z). Together with Eq. (|28|) it also gives a set of three 
equations for Ejj, Ej_j and E^. 
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(2z - S -[■-[■ - E U )(G° - G°_) = GVW (38a) 

G° + + G°_ = GcosA (38b) 

2Z n (G° + - G°_) = -Gsm(£)VW (38c) 

From these equations we obtain G°_ = which implies E- — > oo as well as the three self 

energies. From G°(£' + ) = Gcos(|) and Eq. fl25|) it follows that 

,6. sin 2 (^) / N 

E + = zcosQ + —^ff (39) 
2 Gcos(f) 

We cannot obtain separate expressions for the self energies but GkoV can be calculated 
by knowing the ratios H^/E_, T,^/E_ and E^/EL. From Eq. Q3"8"D and the definition of 
E- (|3"4]), we have, for i?_ — > oo 

(l-tan 2 (^ (40a) 

E n /£L - (cot 2 A - l") (40b) 



E n /£_^ (2cot(|)) (40c) 



Finally the full propagator becomes: 

Z' cot (0/4) -1 ^ 



tan 

r k<7(7' 



2(£+ - e fc ) 



(41) 



y -1 tan(0/4) y 

where E + is given by Eq. (|39|). From Eq. (|29"D we see that the spin-flipped propagators are 
obtained by replacing k by k + Q, i.e. by — e^. 

In Fig. 3 we show the spectral density pk(u)) obtained from these equations for several 
values of (p and k = Q. The Neel case is equal to the local density of states since it is 
k-independent. In particular, no quasiparticle peak is present since the system has no Fermi 
surface. Only the FM background presents a quasiparticle peak that corresponds to the 
dynamics of a free particle. 
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IV. SUMMARY AND CONCLUSIONS 



We have obtained general expressions for the single particle Green's functions describ- 
ing the dynamics of a hole in a generalized spin background in infinite dimensions. Our 
calculations are a generalization of those performed in Ref. 0. 

We characterized the ensemble of spin configurations in a hypercubic lattice by the 
average value of the spin in each sublattice. Closed expressions were found for the symmetric 
case (pAa = Pb-o)i interpolating between Neel order and the unpolarized random case. 

We also performed the calculations considering a magnetic field acting on the system 
that tilts the spins in such a way that spins of different sublattices form an angle <fi between 
them. In this case we obtain expressions for the local and k-dependent propagator as a 
function of the angle 0, so the results interpolate between the FM and Neel configurations. 

The local density of states gets narrower when departing from the Neel order and acquires 
exponential tails that come from the contribution of FM clusters to the hole motion. The 
spectral function p^{uj) also shows this behaviour. A quasiparticle peak is seen only in the 
FM case and it corresponds to the movement of a free particle. 

A self-consistent calculation for finite dimensions can also be carried out in a similar 
way as in Ref. ]I|, considering the exact expression for G° in ci-dimensions and scaling the 
hopping as in Eq. (0). Nevertheless, one has to take into account that in this calculation 
only loop trees have been considered. Paths with loops on which a hole walks around more 
than once or with multiply connected loops, are suppressed with respect to the loop tree 
|T7| by some integer power of 1/d. Also, in the Neel case, at the band edges, where the 
density of states is small, corrections become very large. Another important feature that is 
neglected in these calculations is the presence of spin fluctuations that become important 
at low dimensions. Nevertheless, this approximation can always be corrected by properly 
taking into account all the relevant hole paths (for the 2D Neel case see Ref. fl8f ) or by 
allowing background restoring spin flips along the path. 
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FIGURES 

FIG. 1. Schematic representation of the spin configuration in the tilted background, in the 
presence of a magnetic field h. 

FIG. 2. Density of states for the tilted background for different values of a = cos(0/2): a = 1 
(FM case, dashed-dotted line), a = .8 (dashed), a = .5 (dotted), a = .3 (full), a = (Neel, 
long-dashed line). 

FIG. 3. Spectral density for the tilted background for k = Q and different values of 
a = cos(0/2): The symbols correspond to those of Fig. 2. The FM case presents a delta function 
at lo = 
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Fig. 2 




